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We investigate possible realization of the phantom-like behavior in the framework 
of /(i2)-gravity models where there are no phantom fields in the matter sector of the 
theory. By adopting some observationally reliable ansatz for f{R), we show that it 
^Jj' is possible to realize phantom-like behavior in /(ii)-gravity without introduction of 
phantom fields that suffer from instabilities and violation of the null energy condition. 
Depending on the choice of f{R), the null energy condition is fulfilled in some subspaces 
of each model parameter space. 
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1 Introduction 



One of the most astonishing discoveries of the last two decades is the observation of a 
positively accelerated phase of cosmic expansion. This amazing result comes from several 
sources of observational data such as: measurements of luminosity-distances of supcrnovae 
type la (SNIa) [1], the cosmic microwave background (CMB) temperature anisotropics with 
the Wilkinson Microwave Anisotropy probe (WMAP) sateUite [2] , large scale structure [3] , 
the integrated Sachs- Wolfe effect [4], and the weak lensing [5]. Indeed, general relativity 
with ordinary matter content of the universe leads to a decelerating universe and therefore 
it cannot describe this accelerating expansion which has been confirmed by a huge amounts 
of observational data. In order to realize this late-time acceleration theoretically, several 
approaches have been proposed. One possibility is to consider an extra source of energy- 
momentum with a negative pressure in the matter sector of the Einstein field equations. 
However, the nature of this extra component (the so called dark energy) is yet unknown 
for cosmologists. A very simple and popular candidate for dark energy proposal is the 
cosmological constant [6], but this scenario suffers from some serious problems such as huge 
amount of fine-tuning and coincidence problems. Beside these problems, this scenario has 
not a dynamical behavior because of a constant equation of state parameter (cja = — !)• 
Another suggestion for dark energy is the dynamical models that include various scalar 
fields such as quintessence, k-essence, chaplygin gas, phantom fields, quintom fields and so 
on [7] . On the other hand, one of the most important results of the observational data comes 
from WMAP5 that the equation of state parameter of dark energy can be less than —1 and 
even can have a transient behavior [8] . While general relativity with one scalar field cannot 
realize such a crossing behavior, non-minimal coupling of scalar field and gravity leads to 
this crossing phenomenon [9]. 

There is another approach to realize the cosmic speedup: modifying geometric part of the 
gravitational theory. This proposal can be realized in braneworld scenario ( DGP model and 
its extensions [10]), string inspired scenarios ( Gauss-Bonnet terms in the action [11]) and 
so on. A very popular modified gravity model is the so called /(i?)-gravity [12] where f{R) 
is an arbitrary function of the scalar curvature R. This scenario has the interesting feature 
that choosing an observationally reliable f{R), it is possible to describe the early infiation 
as well as the late time acceleration of the universe in a fascinating manner [12]. Recently it 
has been shown that one can reahze the phantom-hke effect ( increasing of the effective dark 
energy density with cosmic time and an equation of state parameter less than —1) in the 
normal branch of the DGP cosmological solution without introducing any phantom fields 
that violate the null energy condition (NEC) [13,14]. This type of the analysis then has 
been extended by several authors [15]. The main goal in these studies is the realization of 
the phantom-like behavior without introducing any phantom fields in the matter sector of 
the theory. In fact, since phantom fields suffer from instabilities and violate the null energy 
condition, it is desirable to realize this behavior without introduction of phantom fields. 
With this motivation, in this paper we introduce another alternative to realize phantom 
like effect: We study possible realization of this behavior in the framework of /(i?)-gravity 
models. We consider some observationally reliable versions of f{R) gravity and investigate 
the phantom-like behavior of each model without introducing any phantom field that violates 
the null energy condition. Some of these model such as Hu and Sawicki (HS) model have 
passed the solar system tests in a very good manner as well as the perturbation theory [16]. 
We show that all of these models in some subspaces of the model parameter space reahze 
a phantom-like behavior without introducing any phantom fields. We study the conditions 
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that are required in each case to fulfill the null energy condition. 



2 /(i?)-gravity 

In this section we consider the metric formalism of /(i?) -gravity and we summarize the 
field equations of the scenario. The action of a general /(i?)-gravity theory is given by 
[12,17,18,19] 

5 = ^ / d'xV^{f{R) + jCm}, (1) 

where R is the scalar curvature, f{R) is an arbitrary function of R and k = SttG is the 
gravitational constant. The term Cm accounts for the matter content of the universe. Using 
the metric approach, variation of this action with respect to g^i, provides the field equation 

f{R)R^. - \f{R)g,. - V,VJ\R) + g,.uf{R) = /.T^, (2) 

where the prime denotes derivative with respect to R and the matter stress-energy density 
is defined as 

By assuming a spatially flat FRW metric, the Friedmann equation can be written as 



H 



2 



SttG 



Pm 



+ (4) 



where pm is the energy density of the ordinary matter and Pcurv is the energy density of the 
curvature fluid defined as 

Paar. = j7^{^[/(^) ' Rf'^ ' ^HRf'iR)^. (5) 

Throughout this paper we consider the Jordan frame, thus the continuity equation of the 
matter sector can be read as usual 

PM = pM{t = U) = ZHlVtM{l + zf. (6) 

Where VLm is the present day matter density parameter. The continuity equation for the 
curvature fluid is given in the following form [17] 

■ ,_ 3HiQMRf"{R){l + zy 

Pcurv ~r "^-^ \^ ~r ^curvPcurv) [y^J?)]^ ' 

By definition, the pressure of the curvature fluid is given by 

Pcurv = 77^{2^^/"(^) + ^/"(^) + - l[f{R) - RfiR)]]- (8) 
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The equation of state parameter corresponding to the curvature sector of the theory can be 
read as follows 

Rf"{R) + R\Rf"{R) - Hf\R) 



-1 + 



(9) 



\[f{R)-Rf'{R)]-?,HRf"{R) ■ 

From the continuity equations (7) and field equation (4), the Hubble rate can be expressed 
as follows 



2f'{R) 



2>HIQm{1 + zf + Rf"{R) + it Rf"{R) - Hf{R) 



(10) 



where R = 6{H + 2H^). This equation is a very complicated equation and it is very difficult 
to solve it even with the simplest forms of the /(i?) -gravity. 



3 Phantom-like behavior of /(i?)-gravity 

With phantom-like behavior, we mean an effective energy density which is positive and 
grows with time and its equation of state parameter stays less than —1. In this section, 
by adopting some cosmologically viable ansatz, we show that the modified gravity can lead 
to the effective phantom dark energy and phantom-like behavior without need to introduce 
any kind of the phantom (negative energy density) scalar fields that violate the null energy 
condition ( see also [18], [19] for some earlier attempts in this regard). To do this end, the 
modified Priedmann equation (4) can be expressed in a familiar form 
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This relation shows that in f{R) gravity the gravitational constant G can be replaced by 
an effective gravitational constant Gg// = 7^^- "^^^ equation (11) can be recast in the 
following form 



PM + Peff 



12) 



where Pe// — f'{R)Pcurv with p^^rv defined as (5) and therefore effective equation of state 
parameter is given by 



eff 



+ 



Rf"{R) + R 


RfiR) 


- Hf"(R) 




r{R)iWiR)- 


Rf'iR)] 


- 3HRf"{R)] 



(13) 



Now we have all necessary ingredients to discuss phantom-like behavior of /(i?) -gravity. To 
do this end, we consider some observationally reliable ansatz for f{R). 



3.1 Phantom-like effect with f{R) = R + /oi?" 

Phantom-like behavior is the growth of the effective energy density with cosmic time and in 
the same time, the effective equation of state parameter should stay always less than — 1. 
we mean as a first illustrative example, we consider the following ansatz [19] 

/(i?)=i? + /oi?" (14) 
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with constant /o and a. If a < 1, in the small curvature limit the second term dominates. 
Note that in this class of models, a negative a implies the presence of a term inversely 
proportional to R in the action that can lead to the present cosmic speed-up [20]. For a = 
the curvature freezes into a fixed value lead to producing a class of models that accelerate 
in a manner similar to the cosmological constant included models (AC DM). As has been 
stated, for a = —1 this model can describe the acceleration of the universe, but the model 
evolves in the future into an unstable regime where 1 + f'{R) < and it does not contain 
AC DM as a limiting case of parameter space [21]. Now the expression for effective quantities 
in this model take the following form 

„rl SaHRi , , 

Pe// = (l-«)/oi?"[- + -^] (15) 

and 

-1 / R+^(a-2)-HR\ , ^ 

To have an intuition of phantom-like behavior in this case, we adopt the ansatz a{t) — 
aot". It is important to note that this is a solution of the Friedmann equation in our case. 
Especially, for > 1 it gives an accelerating universe which is essentially realizable in f{R)- 
gravity [12]. In table 1 we have shown the acceptable ranges of a to realize phantom-like 
behavior in some subsets of the model parameter space. As mentioned before, theoretically 
negative values of a can account for cosmic acceleration, but they evolve in an unstable 
regime in the future. Especially, it is clear that for a = —1 the null energy condition is 
violated. While for negative values of a, Pe// grows with decreasing z but its values always 
remain negative and the effective equation of state parameter is quintessence-like. The case 
a = is corresponding to an effective cosmological constant with equation of state parameter 
LJeff = —1. Our numerical analysis shows that in this case phantom- like behavior can be 
realized if 0.5 < a < 1. The case a — 1 with a redefinition of the Newtonian gravitational 
constant is corresponding to general relativity. For a > 1, the effective equation of state 
parameter lies in the non-phantom region of parameter space and therefore it cannot account 
for phantom-like behavior. In figure 1 we plot the effective energy density and equation of 
state parameter of the model versus the redshift z ior a — 0.5 (note that this choice is 
corresponding to f{R) = R + fQ\fR model). As figure shows, in this case the effective 
energy density increases with decreasing z and the effective equation of state parameter is 
less than —la typical realization of the phantom-like effect. From figure 2 we can derive the 
acceptable ranges of a to fulfill the null energy condition. As this figure shows, null energy 
condition is respected in some subspaces of the model parameter space and not in the entire 
parameter space, we note that for a < the phantom-like prescription breaks down since 
in this case pe// is negative. 

3.2 Phantom-like effect with \nR term 

In this subsection we consider a modified gravity scenario with In R term in the form [19] 

f{R)^R + f3\n^ + ^R"'. (17) 
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Table 1: Acceptable range of a to have a phantom-like behavior with f{R) = R + /o-R"- 





a < -2 


-2 < a < -1.5 


-1.5 < a < -1 


-1 < a < -0.5 


-0.5 < a < 


a ^ 


< q: < 0.5 


0.5 < a < 1 


a > 1 




negative 


negative 


negative 


negative 


negative 


cosmological constant 


positive 


positive 


positive 




"e.f.f > -1 


UJ^ff > -1 


"e.f.f > -1 


"e.f.f > -1 


^e.f.f > -1 


^e,f,f - -1 


"e.f.f < -1 


"e.f.f < -1 


UJ^ff > -1 




not respected 


respected 


not respected 


respected 


respected 


respected 


not respected 


respected 


respected 




Figure 1: Variation of the effective dark energy density versus the redshift (left hand side). The 
effective dark energy density increases with decreasing z and therefore shows a phantom-like behav- 
ior. The effective equation of state parameter versus redshift (right hand side) which has entered 
in the phantom phase in the past. 




-0.4 - 



-0.5 1 ' ' ' ' ' 1 

-3-2-10123 

a 

Figure 2: pe// ~^ Peff versus a with z = ±1,0. The null energy condition is violated for some 
values of a, but there are subspaces of the model parameter space that respect this condition. 
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The second term in this ansatz containing InR, is growing at small curvature. Basically 
this term is induced by quantum effects in curved spacctimc. It has been shown that this 
model has a well defined Newtonian limit and is able to provide the late time acceleration 
without need to introduce any dark energy component [19]. Choosing m = 2, this model 
leads to a very interesting result: unification of the early time inflation and the late time 
acceleration. On the other hand, considering term can suppress the instabilities arises 
in the perturbation theory of the model as well as improving the solar system bounds, 
consequently the theory can be viable [19]. 

Now, the effective quantities in this model attain the following forms 



Peff 



P 



2/1 



, „9/l 3mHR\ 



SmHR^ , 3pHR 



(18) 



and 



-1 



a;, 



e// 



1 + ^ + ^7^' 



'm—1 



{R - HR) 



+ m{m - l)7i?™"2 + ^ + m{m - l)(m - 2)7/?" 



(1 + ^ + m7i?" 



-2& + fl^# + (l-^)7i?^ 





3mHR\ 


, ZI3HR 




^ ) 





(19) 



In figure 3 we plot the effective energy density and equation of state parameter versus 
the redshift z for m = 2. As this figure shows, in this case the effective energy density 
has a growing behavior with decreasing so it displays the phantom-like behavior without 
introducing any phantom field. The effective equation of state parameter in the late times lies 
in the phantom region with no crossing behavior. To realize phantom divide line crossing we 
can introduce for instance a canonical scalar field in the matter sector of the theory. In figure 
4 we have investigated the acceptable ranges of m to satisfying the null energy condition 
for z = ±0.2, 0. This figure shows that this model respects the null energy condition for 
m > 1.3. As has been pointed out in [19], the presence of higher derivative terms like 
R^ (which may be responsible for early time inflation) in this model helps one to pass the 
existing arguments such as instabilities and solar system tests against such modification of 
the Einstein gravity. 
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1 -0.5 0.5 1 z 1.5 2 2.5 3 ""-^ " 0.5 z 1 1.5 



Figure 3: Variation of the effective dark energy density versus the redshift (left hand side). The 
effective dark energy density increases with decreasing z, so it shows a phantom-like behavior. The 
effective equation of state parameter is less than —1 in the small redshifts ( right hand side). 
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3.3 The Hu-Sawicki model 



One of the most interesting modified gravity model fias been proposed by Hu and Sawicki 
(HS model [22]) that can escape the severe constraint imposed by the solar system tests. 
The form of f{R) in this model is written as follows 



Cl 



f{R)^R-M'' , (20) 

C2[-M2) +1 

where n > and Ci and C2 arc arbitrary dimcnsionlcss constants while Af has the dimension 
of mass. This model yields an effective cosmological constant which generates the late-time 
accelerated expansion [13]. For > M^, equation (20) can be expanded to find 

/(fl)»fi-Mf + M^^l(^)". (21) 

C2 C2 ^ 

In the limit of % —> at fixed ^, this can be reahzed as an effective cosmological constant 

c| C2 ' ° 

Aeff = M'^^ which produces the late time acceleration of the universe. The effective energy 
density in this model is given as follows 



A A 2Ac2(^)" / 6HR\ 3HRAc2n\^ 



^'^f 2 +2C2(^)- + 1^[C2(^)- + 1]A R' J UrJL^n^ll^p2^''W2^ 



ci^r+i'R^ 



R 



(22) 

where A = M^^. The effective equation of state parameter is a lengthy expression and we 
do not write it here explicitly. Figure 5 ( left hand side) shows the behavior of the effective 
energy density versus the redshift for n = 4. Similar to previous cases, the effective energy 
density increases with decreasing z. The effective equation of state parameter is in the 
phantom phase too, but it never crosses the phantom divide line ( figure 5, right hand side). 
Figure 6 shows the acceptable ranges of n to fulfill the null energy condition. Table 2 shows 
the appropriate subspaces of the model parameter space to have phantom-like behavior and 
fulfilling the null energy condition in the HS model. 
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Table 2: Acceptable range of n to have a phantom like behavior. 



value ot 71 


71 < 1.1 


1.1 < 71 < 1.7 


1.7 < n < 2.2 


2.2 < n < 2.7 


2.'7 <n < 3.4 


n > 3.4 




growing but negative 


growing but negative 


decreasing 


decreasing 


growing but negative 


growing and positive 




"e,fi > -1 






^ef f < -1 


"e/ f > -1 


"c.f.f < -1 


null energy condition 


not respected 


respected 


not respected 


respected 


not respected 


respected 



0.5 1 , , , , , , , , , , 1 



0.4 - 




-0.1 0.1 0.2 0.3 0.4 0.5 2 0.6 0.7 0.8 0.9 -0.5 0.5 z 1 15 



Figure 5: The effective dark energy density versus the redshift (ieft hand side). The effective 
dark energy density increases with decreasing z and therefore shows a phantom-like behavior. The 
effective equation of state parameter remains in the phantom region for small values of redshift 
(present day and future times evolution of the universe) (right hand side) . 
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Figure 6: The null energy condition is respected in some subspaces of the model parameter space. 

4 Summary and Conclusion 

In this paper we have studied possible realization of the phantom-like behavior in some viable 
f{R) gravity models. By phantom-like behavior, we mean increasing of the effective energy 
density with cosmic time while the effective equation of state parameter is less than —1. We 
have shown that some models of f{R) gravity can display a phantom-like behavior without 
violating the null energy condition in some subspaces of their model parameter space. To do 
this end, first we have considered a modified gravity model with f{R) = R + /o-R" and we 
found that the phantom-like behavior can be obtained in the region of parameter space with 
0.5 < a < 1 and in this domain null energy condition is respected. Although for negative 
values of a the effective energy density has an increasing behavior with cosmic time, but 
the null energy condition is violated. In the second stage, we have considered a model of 
modified gravity with a Ini? and an additional power law term. With a suitable choice of 
the model parameters, this model which has potential to describe the early time inflation 
and late time acceleration of the universe, accounts for realization of the phantom-like be- 
havior too. There are appropriate subspaces of the model parameter space that null energy 
condition is respected for this choice of f{R)- Finally we have considered the Hu-Sawicki 
which has a very good phenomenology and has been successful to pass the sever constraints 
imposed by solar system tests and the perturbation theory. We showed that this model is 
also capable to account for phantom like behavior without violating the null energy con- 
dition with suitable choice of the model parameters. The main feature of this work is the 
realization of the phantom-like behavior without introducing any phantom fields that violate 
the null energy condition in the spirit of modified gravity. 
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